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ABSTRACT 

This paper deals with the development and analysis of well-posed models and compu- 
tational algorithms for control of a class of partial differential equations that describe the 
motions of thermo- viscoelastic structures. We first present an abstract “state space” frame- 
work and a general well-posedness result that can be applied to a large class of thermo-elastic 
and thermo- viscoelastic models. This state space framework is used in the development of 
a computational scheme to be used in the solution of an LQR control problem. A detailed 
convergence proof is provided for the viscoelastic model and several numerical results are 
presented to illustrate the theory and to analyze problems for which the theory is incomplete. 
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1. Introduction. During the past few years considerable attention has been de- 
voted to the development of smart materials and structures (see [B]). One approach 
to this class of problems is to use shape memory alloys as actuators in active control 
designs. These alloys are best described by thermo- mechanical models consisting of 
coupled (and nonlinear) hyperbolic and parabolic partial differential equations. The 
development of computational algorithms for designing controllers for such systems is 
an immensely complex problem and the subject of several ongoing research projects. 
In addition to the obvious difficulties related to the nonlinearities, the basic thermo- 
elastic coupling often gives rise to nonstandard mathematical models and leads to 
several problems in developing computational algorithms for control. Therefore, the 
computational methods for controlling a linear thermo-elastic system may be viewed 
as a first step toward the ultimate nonlinear problem. With this motivation in mind, 
we consider the problem of controlling a class of coupled partial differential equations 
that describe the linearized motions of a thermo-mechanical structure. The basic 
approach is to combine approximation theory with state space modelling to develop 
convergent computational algorithms for LQR control designs. 

In this paper we consider the questions of well-posedness and convergence of ap- 
proximation schemes for a class of abstract linear systems of the form 

z{t) = Az{t) + Bu(t ), z( 0) = z 0 (IT) 

on a Hilbert space Z. The main concern of this paper is with a general class of partial 
functional differential equations (PFDEs) arising in the modeling of viscoelastic and 



thermo-viscoelastic systems, for example (see [BCLM], [BMC], [MH]), coupled 
equations of the form 

d 2 d T d [° d 1 

a dt? y ( t ' x ^ = Y [Y^ 1 '^ + J_ r 9 ^Y y(<i + s,x ^ ds 

-'y-7^0( t , x ) + b{x)u( t ) i ( 1 . 2 ) 

Wt 0lt ' l} = K §? 0{t ' x) ~ < U > 

where y represents displacement and 6 is the deviation from the reference temperature 
Oq. Equations of viscoelasticity (e.g., (1.2) with 7 = 0) have a special structure which 
has been used by Fabiano and Ito [FI] to formulate a general w T ell-posedness theorem 

and convergence results. Observe that (1.2) with 7 = 0 can be written as 

d 2 3 2 \ f° 1 

a dii y ( t,X ) = Y? Ty{t ^ + J + s,x)ds +b(x)u(t), (1.4) 

or, in abstract form, as 

T V + J g(s)y(t + s)ds =/(<), (1.5) 

where A is a positive definite, self-adjoint, closed linear operator on a Hilbert space 
y. In [FI] Fabiano and Ito consider equations of this form with singular kernels 
(i.e., g 6 L l (— r, 0)) and establish well-posedness when the state space is taken to be 
T>(A l t 2 ) x Y x L 2 (— r, 0; V(A'l 2 )). In this paper, we also consider equations of the 
form (1.5), but the approach we take also applies to the thermo- viscoelastic equations 
(1.2) - (1.3) which cannot be written in the form of (1.5). 

In Section 2 we develop an abstract framework and a generalized well-posedness 
theorem which we apply to the thermo-viscoelastic system (1.2) - (1.3) with zero 
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boundary conditions and to the general viscoelastic system given by equation (1.5). 
Our approach allows a singular kernel, and it also has the advantage that it does not 
require explicit knowledge of the domain of A in order to write down the state 
space. This property can be useful in applications where A 1 ^ 2 is not a differential 
operator. We also remark that our general framework can be applied to certain finite 
delay systems similar in form to the infinite delay systems considered by Miller and 
Desch in [MD], Miller and Desch prove well- posed ness for a class of equations in 
which the kernel is completely monotonic. 

Approximation of such systems generally consists of tvra steps: first approximate 
the spatial variable (e.g., by means of splines) to reduce the system to a hereditary 
differential system on R n , then approximate the “history” or “memory” term (i.e., 
the integral term in (1.5)). In this paper we will use a variation introduced by Fabiano 
and Ito ([FI]) of the averaging scheme considered by Banks and Burns ([BB]) for 
the second stage. The idea of the “AVE” scheme is essentially to approximate the 
kernel g(s) by a step function: partition [— r, 0] into M subintervals and take the 
integral average in each subinterval. Fabiano and Ito show that the approximation 
scheme converges for an L l kernel using a uniform partition of [— r, 0], but they give 
numerical results which indicate that a different partition using a finer mesh near the 
singularity at zero yields much faster convergence. In Section 3 we modify the proof 
given by Fabiano and Ito for singular kernels and a uniform mesh to include singular 
kernels and the non-uniform mesh. 
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Although we prove convergence in this paper only for the abstract viscoelastic 
system, the proof we give can be modified to include the thermo- viscoelastic model. 
In Section 4 we give some numerical results comparing the viscoelastic and thermo- 
viscoelastic models. 

We will use the following notation. For a function g £ Zd(a,6;Z), we denote 
by L 2 (a,b\Z) the set j/ £ L 2 (a,b\Z) f*g(s) ||/(-s)|l 2 ds < ooj. We denote by the 
symbol Hl(a,b\ Z ) the set of all J~P functions which vanish at the left end-point of 
the interval; i.e. , Hl(a,b;Z) = {/ £ H l (a,b\Z) | f(a) = 0}. Similarly, H R (a, 6; Z) = 
{/ £ H l (a , b\ Z) | f(b) = 0}. For a function x : [-r,a) — > X , r, a > 0, the symbol x t 
for t € [0,a) represents the function x t : [ — r, 0] — » X defined by = x(f + s). If A 

is the infinitesimal generator of a Co semigroup T(-) on a Hilbert space Z satisfying 
||T(f)|| z < Me 0t , then we write A 6 Finally, z n -^z means that z n 

converges strongly to z. 

2. Well-Posedness. A standard technique for establishing well-posedness of a 
system governed by a PFDE is to cast the problem in the form of (1.1) and show 
that A generates a C 0 semigroup on Z , for example, by means of the Lumer-Phillips 
theorem (see [P]). We will use the following version of this theorem: 

Theorem 2.1. Let A be a. closed densely defined linear operator on a Hilbert space 
H. If there exists ft € R such that (Ax, x) < ft(x, x) for all x £ T*(A), and ll(\ Q I — A) 
is dense in H for some A 0 > ft, then A is the infinitesimal generator of a C 0 semigroup 
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T(t) on H satisfying ||T’(<)|| < e 0t . 


2.1. A General Theorem on Well-Posedness. Suppose that X , Y, 0 and W 
are Hilbert spaces, and set Z© = X x Y x 0 x W. Let S be a subspace of Y, and 
suppose we have the following linear operators: 


A 0 : V{A 0 ) C Y -+ Y, A, : Z>(A,) C X -* Y, 

Gi : Z>(G]) C 0 — > Y, G 2 : Z>(G 2 ) C Y — 0, 

G 3 : £>(G 3 ) C0-.0, : P(Ci) C W -> Y, 

D : £>(£>) C IV -» IV, i : X —> W, j : 5 -*■ A. 


Define A, G and G by A = AqAj. G = AqGi and G = AqGi, and define Fq by 


V(F e ) = 


e ) € X x 0 x w 

w J 


x € V(Ai),0 £ T>(G\),w £ V(C\), 
A\X G x 0 4* C\W £ 1^(Aq) 


Fq [ 9 
w , 


— A 0 (Aix 4" G\9 A C x iv ). 


Define Aq by 


V(Aq) = 


A& 


( x 

y 

e 

\w 


{ X 

l I ez 9 
V w 

( jy 


y € S n D(G 2 ),0 g D(G 3 ), 1 


w € 2>(0), : 6 V(Fe) 




\ 


Pe 0 
\ w 
G 2 y 4- G 3 0 

\ + Du> / 


Finally, suppose that j is injective and, for A £ p{D) D p(G 3 ), define L\ : 'D(Lx) C 
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X -* Y by 


V(L X ) ={xe 1 Z(j) I (x,\j-'x,(\I - G 3 )- l G 2 Xj- 1 x,(\I - D)-'i\ x) T e V(A e )} , 

L\x = A 2 j~'x - F & \ (A/ - G 3 )~ 1 G 2 Xj~ 1 x ) . 

\ (XI — D)~ l i\x ) 

We are now ready to state the main result of this chapter. 

Theorem 2.2. Suppose 

(1) i and j~ l are continuous, 

(2) T>(A) is dense in X , S C T>(G 2 ) and S is dense in Y , T>(G) fl T>(G 3 ) is dense in 
0, and V(C) D V(D) is dense in W , 

(3) j(S) is dosed in X , 

(4) Fq, G 3 and D are closed, 

(5) for y 6 S, ||C? 2 y||© < k ||iy|| A - for some k > 0, 

(6) there exists (3 € R such that (A q z,z) z& < (3 (z,z) Ze for all z e T>(A@), 

(7) there exists Ao > (3, Ao G p(D) fl p(G 3 ), such that 7 Z(L\ 0 ) is dense in Y, and 

(8) (X 0 I - D)\V(C) n V(D)} is dense in W, and (A 0 7 - G Z )[V(G) n V(G 3 )} is dense 
in 0. 

Then A & is the infinitesimal generator of a C 0 semigroup S e (t ) on Z© satisfying 

IlSeWII ^ e/3 ‘- 

Proof: Set V = V(A)xSx(V(G)r\V(G 3 ))x(V(C)FV(D)). Then V C V(A@) and 

( Xn \ 

T> is dense in Z©, so V(Aq) is dense in Z©. For n = 1 , 2 , . . . , let I 6 V(Aq), and 

\wj 
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A 

i Vn 

suppose ^ 



/Vn\ 

V’n 

In 

\hnj 



as n — ► oo. Then 


y n € S and jy n = <p n — + ip. Since j(S ) is dosed, there exists y € S such that jy - p. 
But j -1 is bounded, so ||y — y|| < ||j -1 || • H V 7 — Sr*™ II + lls/« — J/ll ~ ^ 0 as n 00 ■ 


Therefore, y = y; i.e., y G S and jy — p. Now, ( 0 n ) £ T>(Fq), f 9 n 

\ Wr, ) \W„ 


X 

6 

w 


and F@ I 0 n 

V Wr, 


0 n — ► 0 as n — * oo. Since F& is closed, ( 9 j G V{Fq) and 


F 0 ( 6 ] = 0. Since jy n —* jy and i is continuous, we have ijy n — ► ijy. We also 


w 


have ijy n + Dw n — + h. Thus, Dw n — ♦ h — ijy. But D is closed, so w G F>(D) and 
Dw = h — ijy which implies that ijy + Dw — h. Next, 9 n — + 9, 9 n E V(G 3 ) and 
G 2 y n + G 3 9 n = 7n -» 7- By (5), ||G 2 (y„ - y) || e < k \\j{y n - y)||* = k \\jy n - jy\\ x -* 
0, so G 2 y n — *■ G 2 y which implies that G 3 9 n — » 7 — G 2 y. Since G 3 is closed, 0 G T>(G 3 ) 
and G 3 9 = 7 - G 2 y, or 7 = G 2 y + G 3 0. Therefore, .4© is closed. Finally, let 
(<p\ 

* ±7l(A 0 7-^©);i.e., 

\ h ) 

{p,X 0 x - jy) x + /v’Aoy - f J 


for all 



+ {7, (A 0 7 - G 3 )9 - G 2 y) Q + {h, (A 0 7 - D)w - ijy) w = 0 

/ x \ 


E D(.4e). Let x € T>(L Xo ). Then 


A 0 j~ l x 
(A 0 7 — G 3 ) -1 G 2 A 0 j ! x 
V (A 0 7 - D)~ l iX 0 x J 


= {4>,L Xo x) y = 0 for 


DM©), so ( 0, A„j ! x - F© ( (A 0 7 - G 3 ) 1 G 2 A 0 j ’x 

\ V (A 0 7 - 7})- , a 0 x //y 

all x E P(La 0 ), which implies 0 = 0 by (7). If x = 0, y = 0 and 0 = 0, then 
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< h,(\ 0 I - D)w) w 
X e T>(A), 



= 0 for all w 6 T)[C) fl T>(D), and hence h = 0 by (8). Now for 
€ T>(Aq), so (ip,\ 0 x) x = 0 for all x 6 T>(A). By (2) this implies 


0 / 

that ip = 0 . Finally, for 6 6 V{G)r\V{G^), 



€ T>(Aq), so (7,(A 0 / — G 3 ) 0 ) e = 0 


w 

which implies that 7 = 0 by (S). Therefore, TZ(\ 0 I — Aq) is dense in Z Q , and this 
completes the proof. | 


We wish to apply this theorem to thermo- viscoelastic systems (equations ( 1 . 2 ) - 
( 1 . 3 )) and to the abstract viscoelastic system (equation (1.5)). Before we proceed 
with any examples, however, we will make some general comments about the kernel 
function g and the space W. 


HYPOTHESIS 2 . 3 . The function g satisfies the following conditions. 

( 1 ) 9 € L l {-r, 0 ), 

(2) g < 0 and g' < 0 on [— r, 0), and 

( 3 ) a = t + J° r g{s)ds > 0. 

Set (7a (s) = —-^g(s), and suppose the space X is given. We will take W to be the 
space Lg(—r, 0 ;A') with inner product given by 

(w u w 7 ) w = J g a (s) {w 1 (s),w 2 {s)) x ds. 

Define the operator D by T>(D) = H^— r, 0 ; X), D = The following lemma is 
proved in [FI]. Since it is the crucial step which allows a singular kernel, we reproduce 
its proof here. 
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LEMMA 2.4. The operator D is dissipative in W. 


PROOF: For w E T>(D), 

/ 0 10 \ 1 & 2 

9 o (s) (^ — w(s),w(s)Jx ds = — J ^c(s)q- s IM«)lll ds - 

Let e > 0 and consider 

h = \ J 9o{s)-^ IM*)|& * 

= ^o(-e) IM-e)lli " \j 0 {~r) IK-r)|li ~\J_ r s'M IM 5 )llx ds 

Since u>(-e) = iu(0) - J° t Dw{s)ds = - J° ( Dw{s)ds , by the Cauchy- Schwarz in- 
equality 

M-«)lli < /_° ^ fjM IUM*)lli ds. 

Note that j«(-e) f -~-r = ( 9 °^ 7 ^ ds < e. Thus, we obtain 
J -c 9o\ s ) J—t 

I><\f 9°(s) ||n<»(s)lli 

for all t > 0. Therefore (Dw,w) w = lim^o I t <0. I 
For h E W, if u>(s) = e s f° e~ a h(a)do, then w E V(D) and (I- D)w = h. Since D 
is densely defined, D generates a Co semigroup of contractions on W by the Lumer- 
Phillips Theorem. In particular, D is closed, and A E p(D ) for all A > 0. 

2.2. A Thermo-Viscoelastic System. We consider the system governed by equa- 
tions (1.2) - (1.3) with boundary conditions given by y{t, 0) = y(t, 1) = 0(t, 0) = 
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0(t, 1) = 0. If we set w(t, s, x) = y(t,x) — y(t + s , x) and use g a as defined above, then 
we can rewrite equation (1.2) as 

= y(t,x) + J g Q (s)~-w(t,s,x)ds — -j-0(t,x) + ^b(x)u(t). 

Define the spaces A', Y and 0 as follows: 

X = //q ( 0,1) with {xi,x 2 ) x = a f x[x 

Jo 

Y = L 2 { 0,1) with {yi,y 2 ) Y = a / y x y 2 , 

Jo 

0 = L 2 (0, 1) with (0 u O 2 ) e = ± f 0,0,. 

Vo Jo 

Set Zq = X x Y x 0 x W where IV = L*(—r,0]X) is as above, and let S C Y be 
given by S = //q( 0,1). Define the following operators: 

T>(Ao) = H'(0,l), A 0 y = -y'eY, 

(7 

T>{A l ) = X, A 1 x = x'€Y , 

v{G0 = e, Gi$ = --0 e y, 

a 

V(G 2 ) = H'(0,1), G 2 y = -'yOoy' e Q, 
v(G 3 ) = Hl( o, 1) n H 2 { o, 1), g 3 0 = k0" e 0, 

2>(C,) = H 7 , Cite = j g a (s)w'(s)ds e Y, 

[ia:](s) = x G W, j : //q = $ ~ * X = Hq is the identity operator. 

With the above definitions, we have 

V(A) = H 1 0 nH 2 , V(C) = L 2 g (-r,0-,V(A)), V(G) = H\ 0,1), 
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and the operator .4.© is given by 


V(A e ) = 



y e H%,9 e Hi n // 2 , 

x'-20 + f O _ T g Q ( S )iv'(s)dseH 1 ( ’ 


•4e 


( x\ 


y 


0 


\w) 



( 


f Tx ( x> - 1° + f-r g a {s)v>'{s)ds) 

~l0o y' kO" 


\ 


y + % 


We now verify the conditions of Theorem 2.2. 


( 1 ) Clearly i and are continuous. 

(2) Clearly D(/4) is dense in A r , 5 C X>(G 2 ) and 5 is dense in y, Z?(G) Pi T>(G 3 ) = 

2>(G 3 ) is dense in 0, and T>(C) D T>(D) = 0 ; V{A)) is dense in W. 

(3) Since j(S ) = X , j(S) is closed, 

( 4 ) We already know that D is closed. It is easy to see that G 3 is densely de- 
fined and dissipative, and for any p E 0 there exists 9 6 T)(Gz) such that 
(/ - G 3 )6 = v? (see [K, p. 147]). Thus, G 3 is closed. Let ( ol \ € X>(F©), 


w n 


X 


X r 


X r , 


<9 n ) I <9 , and Fq ( 9 n } = z/> n Observe that F© I 0 n I = 

tt; n / \ w/ \^n / Wn/ 

— (x' n - -9 n + ^( 5 )^( 5 )^). Now, x n -> x in A' implies that < -+ x f 

a ax \ a / 

in y, 0 n ^ 0 in 0 implies that 9 n —> 9 in y, and u; n — ► u> in W implies that 
J! r in 7. Thus, < - l0 n + f° r g a (s)w' n (s)ds 

— > x' — + f° r g a (s)w'(s)ds in F. Since 4o is closed, F© is closed. 

(5) It is easy to check that ||(? 2 !/||© = ~ — ~ llj 2/ II A' ■ ^ = > 
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(6) For z = ( x,y,0,w) T 6 T>(Aq), 

(A e z,z) = a J y'x' + cr J ^ ~ + J 9 c,{s)w'{s)ds^ y 

+ jJ + *0") 0 + J g Q (s ) J y'w' + (Dw, w) 

= k (G 3 0,0) q -f { Dw,w) w < 0 


w 


since G 3 and D are dissipative. 


(7) We will take Aq — 1. Let x £ T){L\). If we set u>(s) = (1 — e s )ix , then 
(/ - D)w — ix, so (I — D)~ l ix = (1 - e s )ix = (1 - e 5 )x. Note that 

^T x ( x/ + ~ G ^~ lx ' + J _ r 9 °^ 1 ~ eS ) xd ^j 


L\x = x 


= x ~-~r (°i x ' + ^—7/ - cy - v 

<j ax V a 


where a,\ = — 
a 


- /° r s(s)(l - e Xs )ds ] = ^ [r + /° r y(s)e As ds] > 0 for A > 
0. Since X C V, we can think of L * as being defined on a subspace of K. 
Observe that V{L } ) = jy € H^( 0, 1) | a,y' + ^(/ - G 3 ) _1 y' € i/ 1 } is dense 
in Y. Define the operators T\ and T 2 as follows: 


V{Ti) = H\ 0,1), r,y = y\ 

V(T 2 ) = Hl( 0,1), T 7 y = y'. 


Note that Xj and X 2 are adjoint to each other. With this notation we can write 
Ly as follows: 


L 1 



0]/ + 


7 2 0q 

a 


(/ - g 3 )- 


1 


T 2 . 
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Since G 3 = G 3 , we have [(/ — G3) 1 ]* = (/ — G 3 ) Thus, 


o. 


i; = / - -t; 


= I--T, 

a 


+ -G 3)- 1 

a 

o,/ + — (/-G3)" 1 

a 


t; 

T 2 = Lu 


that is, is self-adjoint. Now, for y 6 £>(Ti), 

(ii», »)y = {y,y) Y - a j' ± (a,y‘ + ^(/ - G 3 )- v) y 

= llwllv + a /' G,y' + ^(/ - G,)- V) »' 

= Ml + a. toll, + ^ ((/ - GjJ-'y'.y'),, >|l»llv 


since 7 — G 3 > 0 implies that (/ — G 3 ) _1 > 0. Thus, L\ is one-to-one. Hence 
by Theorem 13.11 in [R], Tl{Li) is dense in Y. 

( 8 ) Since V(G?) C V(G), we have (7 — G 3 ) [X>(G) fl P(G 3 )] = (7 — G 3 )£>(G 3 ) = 
0 from ( 4 ). Next, for h £ P(G), if we set tn(s) = e s J° e~°h(a)dcr, then 
w £ 77]j( — r, 0 ; V(A)) = D(C) fl X>(Z>) and (7 - D)w = A. Thus, P(G) C 
(7 — 7)) [£>(G) D V(D)], and V(C ) is dense in W . 

Since ( 1 ) - ( 8 ) hold, Aq generates a G 0 semigroup on Zq. 

2.3. The Abstract Viscoelastic System. Next we wish to apply Theorem 2.2 
to equation (1.5). Since the space 0 corresponds to the heat equation (1.3), we 
first formulate a version of the theorem using only the spaces X , Y and W. Set 
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Z = X x Y x W , and define the operator F by 


ZW = I ( * ) € X x W | x € V(A 1 ),w € P(C,), A,x + C,n> G V(A 0 ) } , 
= Aq(Aix + Ciiv), 


define A by 


V(A) = 


Ml 

V w 


X 

y | €Z 

w 


y e S,w e £>(£>), 


jy 

/(£) I- 

\ijy + Dw 


and define L\ by 


D(Z a ) ={16 7l(j) | (x, Aj-’x, (A/ - J D)- 1 I Ax) r G V(A)} , 
Z. x x = A 2 j~ l x - F ^ A/ _ £)-i;Ax) * 

We have the following special case of Theorem 2.2. 


Theorem 2.5. Assume that all the hypotheses of Theorem 2.2 except those involving 
the space 0 hold. In addition, assume that F is closed, that A satisfies the inequality 
in (6), and that there exists A 0 > 0, A 0 G p(D), such that 1Z(L\ 0 ) is dense in Y. Then 
A generates a Co semigroup S(t) on Z satisfying ||£’(2)|| < e 0t . 


PROOF: Let 0 = {0}, and define G\ — 0, G 2 = 0 (with V{Gi) = Y) and G 3 = 
0. Define V(F Q ) = j ^ 0 j G X x 0 x W | (*) G 2>(F)J , F Q ^ 0 j = f(*), 


and define Aq by T>(Aq) = 


/*\ 
y 
0 

W 


€ Zq 


y J € T>(A) ^ iAqZb = TAT l z e 
w 
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where T : Z -> Z Q is given by T [ y ) = 

w 


y • (and T -1 : Z e -» Z is given by 


j 1 - 1 I ^ I = ( y J )* Note that if we define L\ as before, then T>(L\) — T>(L\) 
w 




(A/- D)~H\x 


= L\x. 


and L x x = A 2 j l x - F & ( 0 i - A 2 ; 

\{\I - D)~H\x J 

It is easy to see that (1) - (S) of Theorem 2.2 are satisfied, so Aq generates a C 0 
semigroup S Q (t) on Z e . Set 5(0* = T~'S e (t)Tz. Then 5(0 ^ a C 0 semigroup on 
Z whose infinitesimal generator is A. I 

We wish to apply Theorem 2.5 to the equation (1.5). If we set te(s) y ys and 


substitute into (1.5) we obtain an equation of the form 


y + q.4 


y + 




= m- 


(2.1) 


We assume that A is a closed, densely defined, positive, self-adjoint, injective linear 
operator. As we shall see in the example considered below, these assumptions are 
easily verified in actual applications. 

A standard technique (see e.g., [FI]) for reformulating (2.1) as an abstract Cauchy 
problem is essentially to set X = V[AF 2 ) where the inner product on X satisfies 
(x!,x 2 ) a - = ^A 1/2 x,, .4 1/2 x 2 ^ and to take the state space to be 

Z = X xY x I 2 (-r,0;A'). 

The approach we take here is similar but it allows more flexibility in the choice of 
state space in that X is not necessarily contained in V, but X will be in one-to-one 
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correspondence with a subspace of Y. We remark that explicit knowledge of the 
square root of A is not required; it is only necessary to know that it exists. 

With this discussion in mind, let S be a subspace of Y containing T>(A), and let 
’? ’) be a symmetric bilinear form on S such that ^(j/ijJ/ 2 ) = (^^Ayi,y 2 ^j whenever 
Vi £ ^{A) and 2/2 € S . Let A be a Hilbert space and j : S — + X a bijective linear 
operator such that : X -» Y is continuous and (x u x 2 ) x = a(j~ x x u j~ A x 2 ). 
Define A : T>(A) C X — > Y by 


V(A) = {x6 X | j~'x € V(A)} , A = -aAj- 1 . 


With W as defined above, set Z = X x Y x W. Then z(i) 6 Z satisfies 


— ^(t) = Az(t) + col(0,/(t),0), 


where A is given by 


V(A) = 


X 

A ( y 

w 


y ) e Z 

w ) 


y € 5, ui € 2>(Z>), 

(* + 9o(s)w(s)ds^ 6 V(A) | ’ 


jy 


A (*+ I-, 9a {s)w(s)ds^ | . 
jy + Dw 


Theorem 2.5. The operator A generates a C 0 semigroup on Z. 
PROOF: Let A, D and j be as above and consider the operators: 

V(C) - W, Cw = J g Q (s)w(s)ds; 

i : X — » W given by [i:r](.s) = x. 


( 2 . 2 ) 
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In order to apply Theorem 2.5, we must show that A can be factored as A — A A- 


Since A : 12(A) C Y — > Y is positive and self-adjoint, it has a positive square root 

A 1 / 2 . Define A 0 and A by 

V(A 0 ) = D(A /2 ), A) = —olA}! 2 
V(A 1 ) = {x£X\ j-'x E V(A 1/2 )}, A = ~A xt2 j~\ 

Clearly, A = AA- Set Ci = A\C . We now verify the conditions of Theorem 2.2. 

(1) Clearly i is continuous. By assumption, j~ l : IZ(j) = X — > Y is continuous. 

(2) Since S 3 12(A) and A is densely defined, S is dense in Y. Suppose x ± 12(A). 
Then for all x E 12(A), 

0 = (x,x) x = cr(j _1 x, j~'x) = a , 

which implies that i^Ay , — 0 for all y E 12(A). But A is self-adjoint and 

one-to-one, so 1Z(A) is dense in Y ([R, Theorem 13.11]). Thus j~ l x — 0 which 
implies that x = 0. Therefore, 12(A) is dense in X . Finally, 12(C) D 12(D) = 
V(AC) fl 12(D) = Hft(—r, 0; D(A)) which is dense in W. 

(3) Obvious since 7 Z(j) = X . 

(4) We already know that D is closed. To show that F is closed, let ) G 

for n = 1.2,..., let (*") -* (*), and f(£) = y. - y as n - oo. Set 
% = )~ l + /°,So(a)t"n(a)ds) and y = j~' (x + f°, 9 i.Mw(s)ds) . Then 

y n E 12(A) and -aAy n = f( Xn ) = y n -* V as n oo. Since A is closed and 



/ 

fO \ 1/2 

7 
• ^ 

VI 

i 

\\*n ~ X\\ X + ( 

J 9 Q (s)ds\ ||u>„ - u;|| w 
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y G P(A) and — a Ay = y. But this implies that x + J_ r g Q (s)w(s)ds G 'D(A) 
and A ^x + /° r g Q (s)u/(s)ds^ = y, so F is closed. 

(5) Does not apply since 7£((7 2 ) Q 0. 

(6) Let | y I 6 D(.4). Then, using the definition of 


w , 


> M 2/ )’( y )/ = {jy,x)x + ( A ( x + [ 9 a{s)w(s)ds'\ ,y 

\ \W J \W ) f 2 \ V J-T / 

+ J 9a(s) (jy + ^w(s),w(s)^ ds 
= J 9ai s ) (^-^ w ( s )iW(s)^ ds < 0 

by Lemma 2.4. 

(7) Again we take A 0 = 1. For x e T>(Li), ( j~ l x j £ V{A), and 

\(l-e 5 )x/ 

L\x = j~ l x — A + J g a {s)(l — e s )xds S j = j _1 x — a^Ax 


where Qi > 0 is as defined above. Define T by 


V{T) = {y G S | jy G Z>(Z,)}, T = I + q x qA. 

It is easy to see that T>(A) = T>(T) and 7Z(T) = 7l(Li). By the Cauchy-Schwarz 
Inequality, for y G V{T), ||Ty|| • ||y|| > | (Ty,y) | = ||y|| 2 + Ql a (Ay,y) > l|j/|| 2 
which implies ||ry|| > |jy||, and so T is one-to-one. Also, T* = T so 7l(T) is 
dense in Y. 

(8) The proof that (/ — D)[D(C) fl T>(D)] is dense in W is the same as the proof 
given in Section 2.1. | 
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EXAMPLE. A viscoelastic shaft with tip-mass. 


Consider a viscoelastic shaft of length / fixed at one end and with a tip-mass at 
the free end. The equation describing the motion of the shaft (see [BMC] where the 
kernel g(s ) is assumed to be in 7/ 1 ) is 

= ^ r^?/(t,x) + y g(s)-^y{t + s,x)ds +b(x)u 1 (t), (2.3) 

while the boundary conditions are given by 


y(t, 0) = 0, (2-4) 

= - + J + s,l)ds +u 2 {t). (2.5) 

Here y is the angular displacement, a is the product of the density of the shaft with 
its polar moment of inertia, r is the product of the shear modulus and the polar 
moment of inertia, I m is the moment of inertia of the tip mass, and the delay r > 0 
is assumed to be finite. Let / = Rx L 2 { 0,/) with 

»^-||(})| I’-'-T’+'jf* 

Then (2.3) - (2.5) can be written as 



19 


Clearly, T>(A) is dense in Y. Let ^ j G 't’(A). Then 

V 5)’ U)) = ~ J‘ ri ’ = /W 

^rol!(;)r 


> 0. 


(5) 2 *(}) - 0) 


Thus, A is positive and 

’ P 2 I m l + aP 

(and A~ l is continuous). Let (^ 2 ) £ Then 


, so A is one-to-one 


J’ )’ ( J 2 j) = $(0lM0 ~ J q 4>"4>2 = 4’Wi 

= MWW) - j o = ((Ji)’^(S))’ 

so >1 is symmetric. Let ^ ~ ^ G K. Define ^(x) = 1 / a4>(Z)d£ + /mil dt. 

Then 6 X>(A), and A = f|Y so 7J(^4) = K. Therefore, by 

Theorem 13.11 in [R], A * = /l. Finally, X>(A -1 ) = 7£(v4) = T, so A~ l is closed, 
and hence A is closed. Now let S = {( j) 6 Y | il> G HIM, 4>(l) = 7 J, and 

define *(($)’($)) = a J o ^ f ° r (Ji)’(S) 65 - Then ’ 5 - ^ and 

, 2 / 2 ) = {aAyi,y^j whenever t/i G X>(^4), y 2 G 5. Let A' = tf[(0,/) with 


is a 


(xi,x 2 ) x = q Jq x\ x’ 2 , and define j : S — * X by j ^ ^ = 4’- Clearly j 

bijective linear operator, and 

(x l ,x 2 ) x = a x\x’ 2 = a a (i" Iac i»J _1 *2) • 

For 1 6 X, ||j -1 x||y = 4nZ 2 (0 + Jq x 2 < II x IIa'- Thus, j" 1 is continuous. 
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Now, define Z = //£((),/) x R. x L 2 (0,l) x L 2 (—r,0]H}^(0,l)) with 


/ 9 
7 

y w 


= a J^iv') 2 + Iml* + v + J 9o{s) dxds » 


and define 


W) 


( <f\ 

7 

W 


V>\ 

7 

0 

VW 

/ 


G Z 


£ A 

<7 3x 


V 


t/>€tfi(0,/), W = 7, 
u? € //^(— r,0; H^OJ)), 

(<p'(x) + J° r g a {s)£w(s,x)ds) G ^(0,/) 

tp \ 


J 


V>'(0 + I-r 9 o(s)£w(s,l)ds 
<P'(x) + f-r 9 a(s)£w(s, x)ds 


Ip + 


dw 

ds 


The operator A generates a Co semigroup on Z by Theorem 2.5. 


3. Approximation. In this section we consider the problem of finding approximate 
solutions to equations of the form (1.1). Let S(t ) be the Co semigroup generated by 
A. We construct a sequence (Z n , P% , A n ) where Z n is a finite dimensional subspace 
of Z, P% is the orthogonal projection of Z onto Z n , and A n generates a C 0 semigroup 
S n (t ) on Z n . We then show that S n {t)P$z -*• S{t)z as n -> oo for all 2 G Z using 
the following version of the Trotter-Kato theorem which follows from Theorem 4.2 in 
[P, Chapter 3]. 

Theorem 3.1. Let A G G{M,f3) be the infinitesimal generator of a C 0 semigroup 
S(t ) on a Hilbert space Z. For n = 1,2, ... , let Z n be a finite dimensional subspace 
of Z such that P n I as n -+ 00 where P n is the orthogonal projection of Z onto 
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Z n . Suppose 

(Hi) A n G G(M,j3) is the infinitesimal generator of a C 0 semigroup S n (t) on Z n for 
n — 1,2, ... , and 

(H2) for all z £ Z , (A I — A n )~ l P n z — ► (A I — A)~ l z a s n — » oo. 

Then for all z £ Z, S n (t)P n z — ► 5(2)2 as ?i — > oo, and the convergence is uniform on 
bounded t-intervals. 


We would like to construct a convergent approximation scheme for the thermo- 
viscoelastic system considered in Section 2.2. Nevertheless, for the approximation 
scheme and convergence proof presented below, we restrict our attention to the ab- 
stract viscoelastic system (equation (1.1)) for which the operator A is given by (2.2). 
The convergence proof for the complete thermo- viscoelastic system involves a modifi- 
cation of the proof we give here and can be found in the thesis [L]. This modification is 
rather technical, and yet is a straight forward extension of the proof we present below. 
Therefore, in order to conserve space we present only the proof for the viscoelastic 
model. If we define the operator A 0 : T>(A 0 ) CXxY—*XxY by 

»Mo)= {(J) € X x y\x£V(A), y e s\, 

+(;)-(&)’ 


then we can write A in the form 
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which suggests a two-stage approximation of A: We first approximate A 0 by discretiz- 
ing the spatial variable, typically by means of spline functions. We then approximate 

— I by discretizing the delay variable. In this paper we will use an averaging scheme 
ds 

for the second stage. We follow the construction given in [FI] except that we do not 
require a uniform partitioning of the interval [ r, 0]. 

Let us now proceed with the first stage of the approximation. Define the bilinear 

form <To(' ? •) ° n A x 5 by 


»■> ((£ )■ (j 2 2 )) = - ”■ 

where a is the bilinear form on 5 discussed in Section 2.3. Observe that for ( ) € 

T>{Ao) and € X x 5, 



Now for each positive integer A r , let X N and Y N be finite dimensional subspaces of X 
and y with y N C 5, and define W N = L)(-r, 0; X N ). We define : X N x Y N - 
X N x Y N by restricting <7o to X N x Y N ; i.e., 

(A 0 V,/), xy = a 0 for ,v N € X N x Y”. 


Now 


set Z N = X N x Y n x IV* and define A N : T>(A N ) C Z N — » Z N by 


V{A N ) = X N x Y n x 

/ 


w I .> 1 = 


-4 2 / 




N 


4 n 


(fi) + L 9 ° ls) ( w ° {s) ) ds 




Jy » + » 


\ 

/ 
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For each positive integer M partition the interval [-r,0] into subintervals [t^ , t^}, 
j = 1 , 2 ,..., A/, where 

= < ■•■ <t? <C = 0. (3.1) 

We will say more later about how the tf are chosen. Set af = tf_ , - f or j = 
1,2,..., M , let Xj { denote the characteristic function of [t™, for j = 2, . . . , M, 

and let x™ denote the characteristic function of [<^,0]. Let i = 0, 1 , , M 

be the usual linear spline functions satisfying B™(tf) = S tj . Define the finite dimen- 
sional subspaces W N ' M and W NM of W by 


W N ' M 




ju; € W 
ju> G W 


M >1 

■w = £ bf* e x N 1 . 


Define D N ' M : W NM -> \V NM by D NM w NM = yu (b^ - bf) x f where 

M 

w x,M _ bf* B and b™ = 0. Define the isomorphism i N ' M : W N ’ M — ► W N,M 
1=1 

by i N ’ M w N ' M = b^Xi*- Now define D N ' M : W N - M ^ W N ’ M by D N ' M = 

D N,M (i N ' M ) . To complete the approximation, set Z N,M = X N x Y N x W N,M , 
and for z N ' M = (x N ,y N , w N ’ M ) T e Z NM , define 


^n,m 2 nm _ 




jy N + £)NM w N,M 
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If W N ' M 


M 

then 


4 N, A/ Z N,M _ 


M 


x K + 

1 — 1 

y N 


a n 

A 0 




M W iW _ U) M 

jy N + E xf 

1=1 «i 


\ 


} 


t M 

where ( g a )<* = 9a{s)ds. 

* 

In order to prove convergence of this approximation scheme, we must impose con- 
ditions upon the spaces A v and Y N and upon the partitions of [— r, 0]. Thus, we 
make the following assumptions: 

(Al) Let P% and Py be the orthogonal projections of X and Y onto X N and Y N , 
respectively. Then P% I x , and P$ Iy where I x and Iy are the identity 

operators on X and Y, respectively. 

(A2) For each positive integer M let n M = {if | j = 0, 1, . . . , M } be a partition of 
[— r, 0] satisfying (3.1), and set A M = {1, 2, . . . , A/}. Then there exist positive 
constants e u e 2 and C independent of M such that \ M = Af U Af where 


A 


M 

1 


{; € a" 


af < rA/‘ (1+Cl)/2 }. 


If j e A**, then (g 0 )™ < jj, and Af contains at most M 1 t2 elements of A M . 
Furthermore, af.j(<7 a )f < {ga) f jof f° r J = 2, 3, ... , M, and if j 6 A 1 , then 


1 , 2 , ... ,j -1 e A? 


M 


* 7* 

EXAMPLE 3.2. Suppose if = — E- for j = 0,1,..., A/. Then af = — for all 
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j E A M , so (A2) is satisfied with = 1, c 2 ,C > 0 arbitrary since = 0 and 

{g a )f<{g a )f. 

EXAMPLE 3.3. Let C = J_ r g a (s)ds, and suppose that the t are chosen so that 

c 

{9ot)Y = — for j = 1,2,..., M . Then it is easy to check that for C\ = e 2 = - the 

J\l l A 

partition IT'^ satisfies (A2) for all positive integers M . 

We will refer to these partitions as the “uniform mesh” and the “non-uniform mesh,” 
respectively. The non-uniform mesh of Example 3.3 is the partition suggested by 
Fabiano and Ito in [FI]. 1 lie convergence proof we give below' is a modification of the 
proof given in [FI]. The major difference is that w'e must handle each estimate in two 
parts: one where the length of the interval is small (in this case our argument is the 
same as that of Fabiano and Ito), and the other where the integral of g is small. 

Let P£' M denote the orthogonal projection of Z onto Z NM . In order to apply 
Theorem 3.1, we must show that P^' M -$-* I z as N,M — » oo. For 2 = ( x,y,w) T E Z , 
P Z ' M z — Py y> Pw My ^j > where P(^’ M is the orthogonal projection of IV onto 

W N ' M . Since we assume (Al) it is sufficient to show that P^' M I w . 

Lemma 3.4. For all h E W, P$ M h -» h as N, M — * OO. 


Proof: For w E IT and t E [— r, 0], set w N (t) = P$w(t), and define w N ' M by 


M 


w NM {t) = £ F °r w € V(D 2 ) 

1=1 


\w — U) 


NM I 


IV < w + w"-*/ uw 


,N 


where the first term on the right-hand side tends to zero by the Dominated Con- 
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vergence Theorem. We write the second term as |ju>^ ||jy ^ ^2 where 

Sj = / _1 9a{s) ||u; N (s) - ^(^i)llx ds for i = ^ 2 . Following Fabiano and 

i€A" ■ 

Ito (but replacing r/A/ , the length of each interval using the uniform mesh, by 
r/A/^ 1+Cl ^ 2 ) we get 


Si < r \\D 2 w\\ 2 w 


1 / t* \ 3 

M • - [ — 7 — 7-7777 ) — ■» 0 as M 
3 \A/( 1+fl )/ 2 / 


00. 


Now, for 5 € [— r, 0 ], 



fO 2 

||u; N (s)|| 2 v < ||u>(s)||* = 

/ Dw(Z)dZ = 

Js X 


i 7"TV V V 7? / 

V< 7 o (0 


vv • 


Thus, 


4<2l / ■" *.w IKWI6 + K(«,)l£] * 

i€A" 

< E fcW * < — * 0 as Af 

9a{- r ) ~t, Jt“ 9c(-r)M 2 


CO. 


Hence, for u> 6 T>(D 2 ), 


d n,m 

w — F H / u> 


IV 


< llu;-UA M ll vv + 




NyM > 


2 jju> -^ N ’ M || VV -4 0 as JV, M — > CO. Since P(D 2 ) is dense in W and 


d NM 

vv 


< 

w 

< 1, 


7 n N,M 7 

/l — *\V ^ 


W 


0 as N, M —* oo for all h G W . 


From Theorem 2.5 we know that A G G( 1 , 0 ). In order to show that A N ' M G C( 1 , 0 ) 
for all N,M, it is sufficient to show that A N ' M is dissipative in Z NM since Z NM i^ 
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where we used the Cauchy-Schwarz Inequality and the inequality 2 ab < a 2 + b 2 , and 
from (A2) the fact that (<7a)f-fi/ a f+i < (9a)f f /&?* for i = 1, 2, . . . , M — 1. Thus, we 
have esta b lis hed the following result. 


LEMMA 3.5. If A N,M is as defined above, then A N ’ M E (7(1,0) for ail N,M. 


All that remains to be done in order to establish convergence is to show that (XI — 
A N ' M )~ l Z — * (XI — A)~ l z for all z E Z. For ReA > 0 and z = f ip ^ , consider 
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the equation (XI — A) 1 


2 = 


X 


, or equivalently, 


y 

w 


>- 

H 

1 

<C2 

II 

(3.2) 

/ fO \ 

(3.3) 

Xy - A ( x + J g a (s)w(s)dsj = ip. 

dw 

Xw - jy - J - = h. 

(3.4) 


From (3.4), w(s) = f° e A < s_c) (jy + h(Z))d£, and from (3.2), jy = Xx - or y - 
Xj~*x — Substituting into (3.3) and using as defined above and the fact that 

\e x ( s ~^d£ = 1 — e A-S we obtain 


A 2 j l x — A^a\x — J g 0 ,(s)(XI — D) 1 [v? — ft(s)] — ?/> + Aj <p. (3.5) 

If we define A(A) = A 2 ; -1 - (r + J° r e xs g(s)ds^ = A 2 ;' 1 - a >. A , then for x € 
T>(A\ 

f 0 

A(A)x = rp + AjTV -A g a (s)(XI - D) l [<p - h(s)]ds. 




Now let Pz’ M z = ( ip N ] , and consider the equation (XI — A N ’ M ) P z ' z 

h N,M 


] N V If we define A N : X N -> Y N by (A N x N ,y N ) Y = -a (j~'x N , y N ), then it 


. w 


y 

N f M 


is easy to see that Aq ^ 


P ( ^ = ( 3 f . r ) for ( X N ) 6 X N x Y N . Thus, we have 


,N 


A' N 


Ax 




,N 
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the equations 


\ N * yv n 
X x -jy =? , 


/ m x 

Ay"-/l''(x' v + £(j t ,)> l ")=^, 
' 1=1 ' 


V' - i/' - K-i ~ Vi 1 ) = A f f for i = 1 , 2 , ... , A/, 


a 


M m 

where w N ' M = £ w^xf* and h N ' M = £ /if'*" From ( 3 .S), 
1 = 1 


M 

E 

t=l 


A + ^K' = i?”" + 7 ^ + ^ 


( 3 . 6 ) 

( 3 . 7 ) 

( 3 . 8 ) 


or 


w. 


M 


(l + q; u A) [ie-Yj + q; w (jy N + A; v/ )] for * = 1 , 2 ,..., M , 


where w ^ 1 = 0 . By induction, w* f = ^ 


i=l 


n(i+<A)-' 


i=k 


°t* 0 y N + h?). From 


( 3 - 6 ), jy N — Ax' N — <,5^ which implies y N = A j 1 x ' v — j 1 < ( ? Ar . Substituting into ( 3 . 7 ) 


we obtain 


A/ 


fc=i L/=* 


tr'^-A" / + »BrfE 

V 1=1 k-l 

M 

= r + a rv w - E 


no+^A)- 

nci+of'A )- 1 




1=1 


By induction, A ^2 

k= 1 


no+^A )- 1 

k~ 1 U=* 

A^ a/ (A) = 


fr=i L/=* 
A/ _ 


“f (v N - A") • (3.9) 


= 1 — 17 (A + oc^f A j 1 . Define 


k = l 


M 


i+B^h-nc 1 +«*'*) 


-i 


i=i 


t=i 


M i fO 

Then, since E(^)f = / s(s)ds = 1, 

i=i a y_ r q 

A"' M (A) = A 1 ;-' - i.A v fr + J° g{s)e M (\, s )ds 


(3.10) 
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where 


•=i \b=i ^ 

Let (A/ - D^) -1 (<p N - /^' M ) = E Cxf- Then 

i=i 

M Mr. 

p N _ = ( A7 _ £ f" Xi y = Y, x C - (Cl - C) X™ 

which implies — hf* = A£f* — ^ ~C) f° r 1 = 1,2, ...,M. Thus, 

l 

(1 + Aq")“‘ ({*_', + a" (V v - /if')], or, by induction, 

tf'-Efnh+w'Hoj' 


kzzl llzzk 


/ 0 M M 

gMZCxt* = Y,(9*)i { C,s o 

r i = l i : = 1 

J° g a (s) (XI - D nm )~' (<? N - h NM ) ds 


*= i l/=* 


AJ l * 

=e^)"E no+^H - ' (3.ii) 


Therefore, from (3.9), (3.10) and (3.11) we obtain 


& N - M (\)x N = 4 > n + X .TV* ~ ^' V /_° $•(*) ( A ' - tD^)" 1 (*>" - A*") (3.12) 

In order to complete the proof of convergence, we must show that x N — ► x, y N — + y 
and -ku as A^, M — + oo. First we need the following lemmas. 


Lemma 3.6. For X > 0, (XI - D NM ) 1 P^' M h — ► (A I - D) l h for all h 6 W. 

PROOF: Let w = (XI - D)~ l h and w N ' M = (XI - D N ' M )~ l P^ M h. Set w NM = 
(i NM )~ 1 w NM , and take w N * M = P^,' M (Dw)d(. It now follows that 

||t N ’ M ( w NM - 53 a ' w )||,v ^ A ||Av M (“> - • (3-13) 
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Now, 


w - i™w NM II, v < |[u> - 


where 


\N,M ||2 _ 


Dw - Pi;r{Dw)\ di 


1 2 

< — Dio - p£' m (Dw) ^ 0 as N , M — > oo 
2 w 


by Lemma 3.4. If we set (Du;) = ]T) CiXi, then we get 


~n,m_ 


M ri M 

= £ / l|f."ll v (* - ( ") 2 * = S. + S 2 . 


To estimate the term Si we again follow Fabiano and Ito (replacing the term r(M 

r 2 

by r/A/^ 1+tl ^ 2 ) to get Si < ||Z?tx;||^, — > 0 as M —* oo. Define the norm |{ - 1[ j 


on W by Hi; = !° r |M«)I& da. Then P™ (D, 


ds = 


M pfL, 


E / "||^|li*=E<*?'inib Forte 6 IV, 


mi; = / M»)iix ds - g ( l_ r ) J 


>) IMMli ds = ^t-y IHIJ, , 


and so, for all j , 


af ikHli < £ K^HL * £<bj • 

Thus, for 5 S ||{f Hi (. - tfY < (of) 2 ||^||’ < jppll 


Hence, 


52 - S ~ 7~ \ H^llvv / 5a (5)^ < 

:rt5a(-r) 


rPHI?yg 

9o(~ r )A/ fj 


0 as M — » 00 , 
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and so ||w - l ' A '- M G) N - M || lv - 0 as JV, M -> oo. Thus, from ( 3 . 13 ) it follows that 
|| W NM _ ;N,M£/v,a/|| iv _» 0 as N, M -» oo. Therefore, 




\w — w 


w — z w 


w 


+ - w N ' M 


w 


0 as iV, M — ► oo. 


f° 

Lemma 3.7. For A > 0, J <7a(s)|e A5 - e A '(A,s)|ds -» 0 as M -» oo. 
PROOF: By definition of e A/ (A,s), 


/ ffa(-s) |e A - s - e M (A,s)| ds = ^ / 5 a(s) 

>/-r i=l A 


- n 0 + «J'A) 

i=! 


-1 


(is 


2C 


= 5 j + S 2 . Let e > 0 , and choose Mq large enough that if M > M 0 , then < c , 


e A(r/A/0 + fl)/2) _ ^ < _ anc J 


e , (rA)V A e 


< For M > Mo, 

2M £l 2 


9C 

M 9 a( s ) ds - < c - 

i6A« 


Let Af = {1,2, . . .,n} and suppose s 6 [< A/ ,0]. Then for some i € Af , tf* < s < 
which implies that e A1 !“ < e« < e A '"., so 0 < e« - e A 'r < e 1 ^ (e A »." - l) 
But, tf* = — J2)=i a f which implies 


< 


e ^_ n(1+< Ar 

j=i 


n^-no+^r 1 


i=i 


j=i 


< 


n (i + a ) _ ri e 

i=i 


\o M 
) 


j= 1 


Now, e Ao J M = (l + Aa^) + -e A ^ (Aa Af ) 2 for some between 0 and a A/ . Set Oj = 
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(l + Aa^) and b 3 — -e Xi > (Aq^) 2 . Then 


ri( i+Q f A )-ii e 




i=i 


i=i 


ek “Iik+m 


i= i i=i 


i t 


j-i 


By induction, FT ~ E[ (F> + &j) = - X>j IT (a* + &*) II «/• Thus, 

j=i i=i i=i t-=j+i /=i 


n(l + <A)-rie-; 


J=1 


i=l 


= E (Aaf) 2 fl ff (> + 


i=i 


i-l 

n 

k-=j+i 1=1 


^ I AfV'M ,M \ 2 / 1 Ar , f ( Ar V _ eAr (-^ r ) 2 C 

- o e ) - 2 e ’ J ' ( ji/(i+e,)/2 ) — 23/ e > < 2’ 

J=1 ' ' ^ 


i-l 


i-i 


since n e x < ft (\ + a?< \) < e x < f\ n t x < = e*? < e x \ Hence, 


t=i+i /=i 


t=i+i /=i 


y<", 

s,< E r *.« 


«6A" '• 

A/ 


„\S \t M 

e — e • 


+ 


e^-ntl + ^A)- 

i=i 




rh-i r 0 

< c Y] / g a {s)ds = tj g a (s)d. 


s. I 


We are now ready to prove the main result of this section. 

THEOREM 3.8. For all z E Z, e* 4 * z — > 5(<)z as N,M — > oo, uniformly on 

bounded t-intervals. 

PROOF: As remarked above, we have only to establish (H2) of Theorem 3.1. Fur- 
thermore, it is sufficient to show that x N — * ;r. If we define the bilinear forms 
//(•,•) on A" by g(x i,x 2 ) = A 2 j~ 1 x 2 ) Y + otx {^i, x 2 ) x 1S easy to verify 

that fi(x x ,x 2 ) = (A(A)a:,,X 2 ) K if X] E F>(A)), and /z M (-,-) on A rA/ by /z Af (xi, i 2 ) = 
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1 1 1 hi III in li 1 1 o 1 m «mp 11 ildm 


(A N - M (A)x 1 ,j _1 x 2 ) r , then for ii,x 2 6 X N , 

\h m (x u x 2 ) - /i(x,,x 2 )| < (y g 0 (s)\e xs -e M (A,5)|^||x 1 || x • ||x 2 || x , (3.14) 

and for x G X N , 

H M (x,x) = A J Hj-'iIIJ, + OA ||x||x • (3-15) 

By (3.5), we have for all u G X, 

/x(x, u) = (ij> + Aj"V,i _1 u) v 

+ a (V 1 J g a {s){\I - D)~ l (s? - h{s))ds,j~ l uj , 
and by (3.12), for all u v c A' v , 

+ a 1 j° g Q (s) {XI - .-)- 1 (*>" - ^’ M (^)) ■ 

Let x N = Pj^x. Then, taking u = u N = x N - x N in the above two equations, we get 

fi M {x N -x N ,x N -x N ) 

< *")| + - *") - rt*,2 N - *")| 

+ Ik - /Ik . Ilr’ll - Ik" - *"ll* + A llr’ll 3 • Ik - /'ll, ■ Ik" - *"ll, 

- (XI - D NM )-'P(t M W -h) is- ||J N - i"lk • 

A 

If we set a¥ = - (t + f° T g{s)e M (X, s)ds\ then by Lemma 3.7, of > 0 for M large 
enough. Estimating the right-hand side of the above equation and using (3.14) and 
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from (3.15) the fact that ax < ft* 1 (x N — x N ,x N — x N ), we get 

I*''-*!, 

- ^ [^ril 2 ^") II** ~ x \\x + (J_ 5a(5)|e A5 -e M (A,s)|) \\x\\ x 

+ Hr 1 II • Ik - + A II;- 1 II 2 • Ik - + ( 0^ ds )' /2 ■ 

■ (XI - D)-\v - h) - (XI - D N ' M )~ l Pw' M (y - A)|| 1 - 0 

II w J 

as TV, M — ► oo by Lemmas 3.6 and 3.7, so \\x — ^ jV || x — > 0 as jV, M — ► oo. | 

We are interested in applying this approximation scheme to the optimal LQR prob- 
lem, but it is well known that convergence of the forward problem (i.e., Theorem 
3.8) is not sufficient to ensure convergence of the gain functionals (see [G], [BIP]). 
Convergence of the adjoint semigroups as well as other properties (uniform stabiliz- 
ability, etc.) play a central role in the development of convergent methods for LQR 
problems. If g(-) £ L 2 ( — r, 0) then Z can be re-normed by using the weight e(s) = 1 
in place of g(s). If Z e denotes the resulting equivalent space, one can establish the 
following adjoint convergence (see [M] for a proof). 

THEOREM 3.9. If g(s) £ L 2 ( — r, 0), then for each z £ Z e , 

et" 4 ’ 1 ( Pz^ z S“(t)z as N,M —+ oo. 

Moreover, this convergence is uniform on bounded t-intervals. 

Recently, I\. Ito has announced [I] a proof of this result that does not require the 
additional assumption that g(-) £ L 2 (— r, 0). In general, the question of preservation 
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of stabilizability uniformly under these approximations is not answered. However, for 
certain special forms of g(*) and for problems with additional damping terms (Kelvin- 
Voigt), one can establish uniform stabilizability. More will be said about this problem 
in the next section. 


4 . Numerical Results. We turn now to an optimal control problem governed by 
the basic thermo-elastic equations with viscoelastic and Kelvin- Voigt damping terms. 


For a rod of length 1 the equations of motion become 


d 2 , d 


9 1 d 

(A + 2/i ) — y(<,x) + - j + s,x)ds - 


dx' 


d 


q(3A + 2n)-^0{t,x) + b(x)u(t) 


(4.1) 

(4.2) 


= X ^ + Qtdx V 

Q3 

where for 0 < 0 < -foo, the term provides Kelvin- Voigt damping. If 

j3 = 0 and g{s ) = 0, then (4.1) - (4.2) become the classical equations of thermo- 
elasticity. If g(s) ± 0 satisifies Hypothesis 2.3, then the integral term provides a type 
of “viscoelastic damping” to the system. 

Gibson, Rosen and Tau [GRT] considered the problem with g{s) = 0 and with 
Dirichlet boundary conditions on displacement and Neumann boundary conditions 
on the temperature. In this case it can be shown that the thermo-elastic model (i.e., 
f3 = 0 and g(s) = 0) has zero as an exponentially stable equilibrium, provided of 
course that one subtracts out the constant temperature distribution. Moreover, the 
elastic and thermal modes decouple, and one can use modal expansions in analyzing 
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this system. Therefore, if /? = 0 and g(s ) = 0, then (4.1) - (4.2) with Dirichlet- 
Neumann boundary conditions is exponentially stabilizable (again, after “subtract- 
ing” the constant temperature states), and the standard LQR control problem has a 
unique solution. The same clearly holds for the corresponding systems with (3 ^ 0 
and g(s ) ^ 0. 

We shall consider the more complex problem governed by (4.1) - (4.2) with Dirichlet 
boundary conditions on displacement and temperature. In particular, we impose the 
boundary conditions 

y{t, 0) = y(t, 1) = 0 = 0(t, 0) = 0{t, 1). (4.3) 

If j3 > 0, then (4.1) - (4.3) is exponentially stable. If 0 = 0 and g{s) satisfies the 
basic assumptions in Section 2 with 

i r° 

e 0 = (A + 2/i) + - j g(s)ds> 0, (4.4) 

then it is known that (4.1) - (4.3) has zero as a globally asymptotically stable equi- 
librium (see [W,pp. 203-210]). It is important to observe that, in general, results of 
this type do not extend to problems in two or more space dimensions. If 0 = 0 then 
it is still not known if (4.1) - (4.3) is exponentially stable. There are positive results 
for infinite delay problems with completely monotone kernels p(-) (see [HW]), 
Recently Hanson [H] has shown that for (3 = 0 and ^(s) = 0, the open loop 
eigenvalues of (4.1) - (4.2) are bounded away from the imaginary axis. However, the 
eigenfunctions for this problem do not form a Riesz basis for the natural state space, 
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and hence the question of exponential stability remains unresolved. 

These theoretical issues can have considerable impact on computational algorithms 
and need to be addressed in order to complete the convergence theory for control 
applications. Some progress has been made for the thermo- viscoelastic problem with 
kernels of the form 


g(s) = —e qa /p\/—s, 


(4.5) 


for q > 0, p > 0 (see [L]), although no results exist for the classical thermo-elastic 
model j3 = 0, g(s) = 0. In summary, if /3 > 0 or g(-) has the form (4.5) with 
(4.4) satisfied, then (4.1) - (4.3) is exponentially stable (hence stabilizable), and if 
/? = 0 = g(s), it is not known if (4.1) - (4.3) is stabilizable. 

We shall present several numerical experiments for LQR control of (4.1) - (4.3) 
with various values of j3 and b when <?(•) is given by (4.5). The numerical constants 
are chosen to be the same as used in [GRT] for an aluminum rod of length 1. In 
particular, 

p = 9.S2 x 10" 2 , A = 2.06 x 10 _1 , 

p = 1.11 x 10 - \ a = 1.29 x 10~ 3 , 

c = 5.40 x 10~ 1 , k = 7.02 x 10“ 7 , 


0 o = 68, q = 30, 

P = 1, r = L 

Observe that for these values, the constant 7 = a(3A+2 p) takes the value 1.0S5x 10 3 , 
and hence the coupling between (2.1) - (2.2) is “small.” Note also that the kernel 
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(4.5) has a weak singularity at s = 0. We selected constants j3 = 1.8 x 10 and 
b = 6.0 x 10 4 so that the first elastic mode has open-loop damping of the same order 
of magnitude as the classical thermo-elastic model. The control input function 6(x) 
is given by 

fl, .4 < x < .435; 

6(x) = 

( 0, elsewhere, 

and is the same as used in [GRT]. 

The LQR problem for the system (4.1) - (4.3) is to choose a control function u Q (t) 
to minimize 

r 00 

J= I [£\t) + u\t))dt (4.6) 

Jo 9 

where £(t) is defined by 

1/2 





' d 

2 r 1 

0 

yl 


dx - f e 0 

Jo 

Tx y(t ' x \ 


dx + 


t fm,xt 

^0 Jo 


dx 


(4.7) 


and y(t,x), 6(t,x) is the solution to (4.1) - (4.3) with initial data 

7/(0, x) = 7/ 0 (x), f t y( 0,x) = i7 O (x),0(O,x) = T 0 (x) 
2 /( 0 , x) - y(s,x) = w Q (s,x). 


(4.8) 


For g(s ) defined by (4.5) and b > 0 such that (4.4) holds, there exists a unique 
optimal control law in feedback form 

u 0 (t) = t 0 J -^K l {x)-^y{t,x)dx + p J K 7 {x)J- t y{t,x)dx-^Y Q J q K * {x)0{t,x)dx 

°-kir 


***<*•*> 


0 0 

—y(s,x) - — y(t + s,x)\ g{s)dxds. (4.9) 
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If g(s) = 0 and /? > 0, then again there exists a unique optimal control law in 
feedback form (4.9) with g(s) = 0 = /\ 4 (s,x). If 0 = 0 and $(s) = 0, then the 
most one can say is that if an optimal control exists, then it will have the form (4.9) 
with g(s) = 0 — I\ 4 (s,x ). The existence of the optimal control has not yet been 
established. We had hoped that the numerical results presented below would shed 
some light on this question. However, as will become evident from these numerical 
experiments, we seem to have raised more questions than we have answered. 

The objective of the computational scheme is to produce numerical approximations 
of the optimal functional gains A'i(*), A^(’)> A 3 ( ' ) and A^*,*). In particular, the 
general idea is to introduce an approximation method (such as in Section 3), and 
then to use this method to compute approximating (sub-optimal) gains (■), K 2 (*)> 
K?(') and A^(*,-), The basic questions are i) Do the optimal gains exist? and n) 
Does the particular approximation scheme lead to convergence of the sub-optimal 
gains; i.e., does A"/ v (-) — > A '*(•), i = 1,2, 3, 4 as N — * 00 ? 

The computational results presented below are based on the numerical approxima- 
tion scheme presented in Section 3. We used the non-uniform mesh as described by 
Example 3.3 for the thermo- viscoelastic model. In particular, all runs presented below 
were based on M = 8 subdivisions of [—1,0], and the corresponding approximating 
finite dimensional system becomes 

z N {t) = A N '*z N (t) + B N * s u(t) (4.10) 
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where B N ' 8 = P^' S B and B : R 1 — > Z is defined by = (0, b(x)u, 0, 0) T . The results 
are the same for M = 16 sub-divisions of [—1,0] (a nice feature of the non-uniform 
mesh algorithm). 

In the cases where i) holds, Theorems 3.8 and 3.9 imply convergence of the func- 
tional gains provided one can establish that this scheme preserves stabilizability uni- 
formly in N (see [G]). For the thermo- viscoelastic model with </(•) given by (4.5) and 
the thermo-elastic model with /? > 0, preservation of stabilizability can be shown. 

In order to see the effect of the damping models on the open-loop system, in Figure 
1 we plot the open-loop poles for the first 8 elastic modes using only thermal damping 
(i.e., j3 — 0, g(s) = 0), thermo-viscoelastic damping (/? = 0, g(s) given by (4.5) and 
6 = 6.0 x 10 4 ) and Kelvin- Voigt damping (/? = 1.8 x 10”'), respectively. Observe 
that the damping is extremely small in all three cases. On the other hand, Table 1 
shows that at the higher frequencies the damping provided by the thermo-viscoelastic 
model is an order of magnitude more than the thermo-elastic model, and the damping 
predicted by the Kelvin- Voigt model is two orders of magnitude greater than the 
thermo-elastic model. 
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Elastic 

Mode 

Thermo-Elastic 
0 = 0 = g(s) 

Thermo- Viscoelastic 
6 = 6. Ox 10 4 

Kelvin-Voigt 
0 = 1.8 x 10" 7 

1 

-0.0000054 + 6.5Sz 

-0.0000091 +6.5S?' 

-0.0000093 + 6.58? 

2 

— 0.0000060 + 1 3 . 1 7 z 

— 0.00001S4 + 13.17? 

-0.0000216 + 13.17?' 

3 

— 0.0000071 + 19.79? 

-0.0000325 + 19.79? 

-0.0000423 + 19.79? 

4 

-0.0000086 + 26.47z 

-0.0000501 + 26.47i 

-0.0000715 + 26.47?' 

5 

—0.0000106 + 33.20? 

-0.0000699 + 33.20? 

-0.0001095 + 33.20? 

6 

—0.0000129 + 40.02?: 

-0.0000905 + 40.02?' 

-0.0001566 + 40.02? 

7 

-0.0000157 + 46.93?' 

-0.0001113 + 46.93? 

-0.0002133 + 46.93? 

8 

— 0.00001S9 + 53.95? 

-0.0001317 + 53.95?' 

-0.0002800 + 53.95?' 

9 

—0.0000224 + 61.10? 

-0.0001514 + 61.10? 

-0.0003574 + 61.10?' 

10 

-0.0000263 + 68.39? 

-0.0001705 + 68.39? 

-0.0004460 + 68.39? 

15 

-0.0000491 + 107.5? 

-0.0002565 + 107.5? 

-0.0010S66 + 107.5?' 

20 

-0.0000662 + 151.5? 

-0.0003257+ 151.5? 

-0.0021294 + 151.5?' 


TABLE 1. OPEN LOOP POLES; SMALL DAMPING 

We also considered cases where b — 6.0 x 10 2 , b = 6.0 x 10 1 and 0 = 1.8 x 10 5 , 
0 = 1.8 x 10 -4 . The models with 0 = 0 and b = 6.0 x 10 2 and 0 = 1.8 x 10 4 , 
g(s) = 0 show considerable increases in open-loop damping (especially at the higher 
modes) as illustrated in Table 2. 
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Elastic 

Thermo- Elastic 

Thermo- Viscoelastic 

Kelvin-Voigt 


Mode 

II 

0 

II 

6 = 6.0 x 10 2 

0 = 1.8 x 10" 5 

-* 

1 

-0.0000054 + 6.58* 

-0.0003812 + 6.57* 

-0.0003934 + 6.58* 

f 

2 

—0.0000060 + 13.17* 

-0.0012388+ 13. 16* 

-0.0015617+ 13.17*' 

l 

3 

— 0 .0000071 + 1 9.79i 

-0.0025453 + 19.78*' 

-0.0035214 + 19.79*' 

I 

4 

-0.0000086 + 26.47* 

-0.0041601 + 26.45*’ 

-0.0062913 + 26.47*' 

| 

5 

—0.0000106 + 33.20* 

-0.0059423 + 33.19* 

-0.0098982 + 33.20*' 


6 

—0.0000129 + 40.02* 

-0.0077757 + 40.00*' 

-0.0143768 + 40.02*' 

- 

7 

—0.0000157 + 46.937 

-0.0095789 + 46.91* 

-0.0197699 + 46.93*' 

- 

8 

—0.0000189 + 53.95a- 

-0.0113034 + 53.93*' 

-0.0261290 + 53.95*' 

= 

9 

—0.0000224 + 61.10* 

-0.0129269 + 61.08* 

-0.0335142 + 61.10*' 

- 

10 

-0.0000263 + 68.39* 

-0.0144449 + 68.37*' 

-0.0419944 + 68.39* 

_ 

15 

-0.0000491 + 107.5* 

-0.0207916 + 107.5*' 

-0.1038004+ 107.5*' 

x 

20 

-0.0000662 + 151.5*' 

-0.0260228+ 151.5*' 

-0.2063817+ 151.5*' 

- 


TABLE 2. OPEN LOOP POLES; MEDIUM DAMPING 



The LQR problem for the approximating system (4.10) is solved (by Potter’s 
method) and approximating functional gains K^(-), K^(-), (•) and are 

constructed by the standard Galerkin scheme (see [M] for details). Convergence of 
these sub-optimal functional gains as N —* 00 to the optimal gains can be established 
for 0 > 0 and g(s) — 0 by methods similar to those in [G]. If g(s) is defined by (4.5), 
then this convergence can be established by a modification of Ito’s recent results [I]. 
Nothing is known about the classical thermo-elastic problem. We treat this problem 
numerically. 

Example 4.1. In this problem we set g(s) = 0 = 0 and compute the gains 
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!<?(■) and A'^(-). Note that the finite dimensional system (4.10) is exponentially 
stable (all open-loop poles are in the left half plane) so that these sub-optimal func- 
tional gains do exist. Figures 2, 3 and 4 contain the plots of A 0 (■), i = 1,2,3 for 
N = 24, 28, 32 and 36. Although these plots are similar to those found in [GRT], 
and for the boundary conditions in [GRT] one can prove convergence of the gains, 
it is not clear from Figures 3-4 that these gains are converging. We will return to 
this issue later. 

Example 4.2. In this problem we consider the thermo-viscoelastic problem defined 
by (3 = 0 and b = 6.0 x 10 1 , b = 6.0 x 10 2 . Figures 5, 6 and 7 contain the plots of 
A'N(-), i = 1,2,3 for N = 24, 28, 32 and 36 where b = 6.0 x 10 1 . Figures 8, 9 and 
10 contain the same plots for b = 6.0 x 10 2 . Observe three important features: t) 
For b = 6.0 x 10 1 the system is more heavily damped than for b = 6.0 x 1 0 2 ; ii) the 
functional gains for the case b — 6.0 x 10 1 are smoother than for b — 6.0 x 10 2 ; and 
iitj the convergence of the functional gains is faster in the problem with the most 
open-loop damping. The plots of (•, •) for b = 6.0 x 10 1 and N = 24, 28, 32 and 
36 are shown in Figures 11 - 14. Note that these gains converge rather rapidly. This 
was typical of all the thermo- viscoelastic runs. 

Example 4.3. In this problem we consider only the thermo-elastic model with 
Kelvin-Voigt damping so that (?(s) = 0 and = 1.8 x 10 -4 and 0 = 1.8 x 10 5 . 
Figures 15, 16 and 17 contain the plots of K0 ( • ) , i = 1,2,3 for N = 24, 28, 32 and 36 
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for (3 = 1.8 x 10“ 4 , and Figures 18, 19 and 20 show the same plots for = 1.8 x 10“ 5 . 
Observe that as in Example 4.2 above, the damping affects the smoothness and rate of 
convergence of the functional gains. Also note that when (3 = 1.8 x 10“ 5 the first elas- 
tic mode has the same damping factor (i.e., open-loop pole is — 0.00039 -f 6. 58i) as the 
problem in Example 4.2 with b = 6.0 x 10 2 (i.e., open-loop pole is —0.00038 + 6.5Sz). 

Figures 21, 22 and 23 contain the plots of A' 32 (-) for the thermo-elastic and thermo- 
viscoelastic models and the thermo- elastic model with Kelvin- Voigt damping for var- 
ious values of b and /?. The functional gains for the thermo- viscoelastic problem 
“converge” to the functional gains for the thermo-elastic problem as b — ► oo. The 
same convergence applies as the Kelvin- Voigt parameter /? — ► 0. As illustrated in 
Figure 23, the functional gain A" 32 (’) computed by using the thermo-elastic model is 
the same as the gain computed by using the same model with Kelvin-Voigt damping 
parameter f3 = 1.8 x 10“ Likewise, if 6 = 6.0 x 10 4 in the thermo- viscoelastic model, 
then the functional gain AT(-) is also identical to the gain computed by using the 
thermo-elastic model. This “basic” observation applies to A'? 2 (-) for i = 1,2, 3 and for 
Kf 4 (-) also. In particular, as shown in Figure 24, the functional gain A'f 4 (-) computed 
from the thermo-elastic model is remarkably close to the functional gain computed by 
adding Kelvin-Voigt damping. Although convergence of the functional gains for the 
thermo-elastic model with Kelvin-Voigt damping is assured by theory, Figures 23 and 
24 illustrate that this convergence can be extremely slow if the damping parameter 
l 3 is small. 
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Figures 23 and 24 also raise the possibility that even in the thermo-elastic problem 
(no structural damping) where there is no proof of convergence, the gains might 
converge if N is sufficiently large. The jagged nature of occurs because of the 

small damping. Moreover, comparing the open-loop and closed-loop poles, one sees 
that the optimal feedback law introduces considerable damping at all frequencies. 

Figure 25 compares the first 8 open-loop poles (+) to the closed-loop poles (*) for 
the thermo-elastic model. It is interesting to note that the same pattern holds for 
all of the higher poles (see Table 3 below), and the thermal damping alone does not 
appear to be a major aid in controlling the higher modes. 


Mode 

Open-Loop 

Closed-Loop 

1 

-0.0000054 + 6.58? 

-0.1078901 + 6.58?' 

2 

-0.0000060 + 13.17? 

-0.0552142+ 13.17? 

3 

-0.0000071 + 19.79? 

-0.0792251 + 19.79?' 

4 

-0.0000086 + 26.47? 

-0.0953868 + 26.47? 

5 

-0.0000106 + 33.20? 

-0.0298920 + 33.20? 

6 

-0.0000129 + 40.02?' 

-0.1097076 + 40.02?' 

7 

-0.0000157 + 46.93? 

-0.0264249 + 46.93?' 

8 

-0.0000189 + 53.95?' 

-0.0946583 + 53.95? 

9 

-0.0000224 + 61.10?' 

-0.0742450 + 61.10?' 

10 

-0.0000263 + 68.39? 

— 0.0549S55 + 68.39?' 

15 

-0.0000491 + 107.5? 

-0.0707290 + 107.5? 

20 

-0.0000662 + 151.5? 

-0.0731157 + 151.5?' 

TABLE 3. 

OPEN-LOOP VS. CLOSED-LOOP POLES 


THERMO-ELASTIC MODEL 
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As illustrated by Figures 26 and 27, the first 8 closed-loop poles for the moderately 
damped (Kelvin- Voigt) thermo-elastic and the thermo- viscoelastic models are essen- 
tially the same as the corresponding closed-loop poles for the thermo-elastic model. 

In view of these numerical results, we first conjectured that although no theoretical 
results exist to prove the existence of the optimal feedback gain for the classical 
thermo-elastic model (with boundary conditions (4.3)), existence and convergence of 
the suboptimal gains do hold. Since the damping for the realistic model is so small, 
the numerical evidence provided by the previous plots is not strong. In order to test 
this conjecture we investigated the convergence of the gains for the problem with 
“artificial” parameters. In particular, we considered the equations in dimensionalless 
form 

d 2 d 2 d 

dp y ^ t,x ' ) = ~ + 

with 6(1) as above and c 2 = 7 = 1.0. The uncontrolled system has considerable damp- 
ing. The first two elastic modes have eigenvalues —.0939 -f 3.15i and —.4066 + 6.47z, 
respectively. The functional gains for this set of parameters are slightly smoother 
than the corresponding gains for the aluminum rod. Figures 28, 29 and 30 show the 
“convergence” of these gains for N = 24,28,32 and 36. It appears that the scheme 
may produce convergent gains (at least for and {■))] however, the behaviour 

of K^i') ' s not clear-cut. Based on our numerical experience, we do conjecture that 
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the classical thermo-elastic model with Dirichlet boundary conditions is in fact stabi- 
lizable. However, we also conjecture that the computational scheme does not produce 
convergent gains (we believe that this scheme does not preserve stabilizability uni- 
formly). We have not been able to provide proofs of these conjectures, and we leave 

them as open problems. 
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K2, Kclvin-Voigl Damping, bcla = 1 .8c-5. N = 24, 28, 32, and 36 
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Figure 25 


Open vs. Closed Loop. N = 32, Thcrmo-Viscoclastic, b = 600 
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Open vs. Closed Loop, N = 32, Kclvin-Voigt Damping, beta = 1.8c-5 



Figure 27 


xlO -3 Kl, Gamma=1.0, N = 24, 28, 32 and 36 
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